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The aim of this paper is to describe the self-propulsion of a micro-robot (or micro- 
swimmer) consisting of N spheres moving along a fixed line. The spheres are linked to 
each other by arms with the lengths changing periodically. For the derivation, we use 
the asymptotic procedure containing the two-timing method and a distinguished limit. 
We show that in the main approximation, the self-propulsion velocity appears as a linear 
combination of velocities of all possible triplets of spheres. Velocities and efficiencies of 
three-, four-, and five-swimmers are calculated. 

The paper is devoted to H.K.Moffatt, who initiated the author's interests in low- 
Reynolds-number fluid dynamics. 
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1. Introduction and formulation of problem 

1.1. Introduction 

The studies of micro- robots represent a flourishing modern research topic, which creates a 



fundamental base for modern applications in medicine and technology, see Purcell (1977) 
Becker et.al. (2003)ljNajafi fc Golestanian (2004)||Dreyfus et.al (2005)[|Earl et.al. (2007) 



Chang et.al. (2007)[[Alouges et.al. (2008)[|Gilbert at.al. (2010)[|Golestanian fc Ajdari (2008)[ 



Golestanian fc Ajdari (2009)||Alexander et.al. (2009)[|Leoni et.al. (2009)[|Lauga (26TT) 
At the same time, the simplicity of both time-dependence and geometry represents the 
major advantage in the studies of micro-robots (in contrast with extreme complexity of 
self-swimming microorganisms, e.g. Pedley & Kessler (1987) Vladimirov et.al. (2004) 



1. |Pedley fc Kessler (1987) ^ 
Pedley (2009)| Polin et.al. (2009)); which makes it possible to describe the motions of 
micro-robots in greater depth 



In this paper, we generalize the theory of a three-sphere micro-r obot of |Najafi fc Golestanian (2004)[ 
Golestanian & Ajdari (2008) to an A^-sphere micro-robot. We employ two-timing method 



and distinguished limit arguments, which lead to a simple and rigorous analytical pro- 
cedure. Our calculation of the self-propulsion velocity of an N-sphere robot shows that 
it represents (in the main approximation) a linear combination of velocities due to all 
possible triplets of spheres. The velocities and LighthhTs swimming efficiencies of three-, 
four-, and five-sphere robots are calculated as examples. This paper represents further 



studies of a problem introduced in Vladimirov (2012a) 



1.2. Formulation of problem 

We consider a micro-robot consisting of N rigid spheres of radii R*, i = 1, 2, . . . JV with 
their centers at the points x*(t*) of x*-&xis (x* +1 > x*), t* is time, asterisks mark 
dimensional variables and parameters. The spheres are connected by N — 1 arms/rods, 
such that the distances between the centres of neighbouring spheres are Z* = |a;* +1 — x* |, 
a = 1, 2, . . . N — 1, see FigQ] (in this paper latin subscripts change from to N, while 
greek subscripts change from to A — 1). We accept Stokes's approximation where masses 
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Figure 1. N = 4 spheres, linked by arms of periodically changing lengths. 



of spheres and arms are zero; the arms are so thin (much thinner than any R*) that their 
interactions with a fluid are negligible. The equations of motion can be written as 



F* 



N 



F*=0 

* r • * 

- — Kj \Xj 

— T* — I* 

1 x a _ l Q 



J}, 



b ik 



(1.1) 

(1.2) 
(1.3) 
(1.4) 



where = 6nr]R*, 7/ is viscosity, dots above the functions stand for d/dt*, and summa- 
tion convention over repeating subscripts is not in use. The forces /* are exerted by the 
arms on the i-th sphere, while F* represent viscous friction. In order to derive (|1.2[) we 
use a classical expression for Stokes's friction force as well as the explicit expression for 
fluid velocity for a moving along x*-axis sphere. The x*-component u* of this velocity at 
distance r* along x*-axis is 

u* ~ 3R*U*/(2r*) 

where R* and U* are the radius and velocity of a sphere, see Lamb (1932) Landau fc Lifshitz (1959)] 
Moffatt (1996) Equality (|1.4j) follows from the fact that the external forces exerted on 
each arm are negligible. The geometrical configuration of a micro-robot is determined by 
given functions 



l * a = L* a + l* a (T) 



(1.5) 



where L* a are mean values and ^*(r) are oscillations, which represent 27r-periodic func- 
tions of a dimensionlcss variable r = uj*t*; uj* is a constant. Since all /* are given, 
then conditions (jl.3|) can be considered as geometrical constraints. Equalities ()l.ll) - (|1.5|) 
represent the system of 2N equations for 2N unknown functions: 

x*(n = (xun, xun ■ ■ /*(**) = (f*(n mn mi) 

These equations contain three characteristic lengths: distance L* between the neighbour- 
ing spheres, radius R* of spheres, and amplitude A* of arm oscillations; the characteristic 
time-scale is T*: 



N 



N 



a=l 



N 

E 

a=l 



max^/iV, T* = 1/lu* 



The dimension of n* can be eliminated from the equations by division of (|1.1[) by k*, 
hence it does not play any role in scaling. We choose dimensionless variables and small 
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parameters as: 

x* = L*x, R* = R*Ri, T* a = X*T a , t* = T% f* = &Trr}R*L*fi/T*- : 
e = X*/L* «1, S ee 3R*/(2L*) «C 1 

Then the dimensionlcss form of f|l . (|1 .5[) is 

RiXit — S 2^ RikXkt/hk — fi, hk = ijfc + skk (1-6) 
kjii 

J2fi = f-I = 0, !,...,!) (1.8) 



where Rik = RiRk, subscript t stands for d/dt, 



k-l 



lik — ^2 In, f° r k>i + l, with Iki = hk otherwise, (1-9) 

n—i 

and similar definitions for La, and 1^. (for example, L13 = L\ + Li etc.). Eqn. (jl.6p can 
be rewritten as (N x iV)-matrix form: 

N 

Ax t = f or ^A 4fc x fct =/. t (1.10) 
fe=i 

a a I Ri for i = fc, /1 1 1 \ 

A U { SR ik I, for^fc 



2. Two-timing method and asymptotic procedure 

2.1. Functions and notations 
The following dimensionless notations and definitions are in use: 

(i) s and r denote slow time and fast time; subscripts t and s stand for related partial 
derivatives. 

(ii) A dimensionless function, say G = G(s,r), belongs to class T if G = 0(1) and 
all partial s-, and T-derivatives of G (required for our consideration) are also 0(1). In 
this paper all functions belong to class I, while all small parameters appear as explicit 
multipliers. 

(iii) We consider only periodic in t functions {G 6 V : G(s, r) = G(s, t + 27r)}, where 
s-depcndcncc is not specified. Hence all considered below functions belong to V{\1. 

(iv) For arbitrary G S V the averaging operation is: 

(G} = ±- £° +27T G( s ,T)dr = G( S ), Vr (2.1) 

(v) The tilde-functions (or purely oscillating functions) represent a special case of V- 
functions with zero average {G } =0. The bar-functions (or mean-functions) G = G(s) 
do not depend on r. A unique decomposition G = G + G is valid. 
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2.2. Asymptotic procedure 
The e-dependence of kk f| 1 -6[) leads to the presentation of matrix A (jl.lip as a series in e 

A = B + eSA' Q + ..., B = A + 5Ai (2.2) 



= diag{i?i,i? 2 , ...,Rn}, 



for i = k, 

Rikhk/L^k for * k 



where we do not present the expression for Ai since it does not affect the final answer. 
In double series (with small parameters e and S) matrices Ao, Ai and A appear in the 
terms of orders £°S°, e°5 1 , and e 6 . 

The introduction of fast time variable r and slow time variable s represents the crucial 
step of our asymptotic procedure. We choose: 

t = t, s = e 2 t 



This choice can be justified by the same distinguished limit arguments as in Vladimirov (2012b) 



Here we present this choice without proof, however the most important part of this proof 
(that this choice leads to a valid asymptotic procedure) is exposed and exploited below. 
After the use of the chain rule 

d/dt = d/dr + e 2 d/ds (2.3) 

we accept (temporarily) that r and s represent two independent variables. The substi- 
tution of (|2.3[) . (|2.2[) into (| 1 . 10[) gives its two-timing form: 

(I +^A + ...)(a; T +£ 2 ;r s ) = / (2.4) 

Unknown functions are taken as regular scries 

x(t,s) =x (s)+ex 1 (T,s) + f(T,s) = f (r,s) + ef 1 (T,s) + ... (2.5) 

where Xq = 0, which means that long distances of self-swimming are caused by small 
oscillations: |x | ^> \ex\(t,s)\. 

2.3. Successive approximations 

The successive approximations of (|2.4j) . (|2.5|) yield: 
Terms of order e = 1: / = 0; 

Terms of order e : BoCCo T = fx, its average gives f 1 = and the oscillating part is 

I Sir=7i; (2-6) 
Terms of order e 2 : B0S2T + SA' Xi T + Boccbs = / 2 ; its averaged part is 

B ^ 0s + S(A' x lT ) = / 2 (2.7) 
Force f 2 can be excluded from (|2.7[) . (|1.8[) as: 

I-m x 0s + 5I- (A 5i T ) = (2.8) 

The averaged self-propulsion motion means that the rate of changing is xois = Xq s with 
the same function Xq(s) for all spheres; therefore we write Xq s = Xq s I. Hence (|2.8[) gives 

X 0S ^ S 1 -^ (2.9) 
I ■ AqI 

where in the denominator matrix Bo is replaced with Ao, since we consider only the main 
(linear in S) term in (|2.9|) . Expression (|2.9j) still contains unknown functions x\ T which 
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can be determined from (|2.6j) with the use of constraints f|l .7)1 . f)l .8[) . Indeed, equation 
(|2.6|) (with the terms required for linear in 5 precision in (|2.9j) ) gives: 



xi r = (A ) 1 f 1 =g, g = (g 1 ,g 2 ,...,g N ) 

gi = fu/Ri, (So)- 1 = diag{l/i?i, 1/R 2 , 1/Rn} 



(2.10) 



One can see that (|2. 101) . (|1 .7[) yield x a +i, T — x a ^ T = g a +i — g a = lar, while (| 1 . 8|) leads to 
Si R*9i = 0- Both restrictions can be written with the use of N x N constraint matrix 
C: 



Cg = l T , 
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The substitution of its inverse form 



It 



(2.11) 



into ([2J]) yields 



where 
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Pn, 
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(2.12) 



(2.13) 



The expression for inverse matrix (|2.13j) can be checked by direct calculations of product 
C _1 C; these calculations become particularly simple if the matrix in the right hand side 
of (|2.13[) is decomposed into two matrices: one containing all p's, and another with iden- 
tical rows (pi, p2, P3, ■ ■ ■ ,pn-i, 1)- It is worth to emphasise that the presented analytic 
procedure is especially simple, since in order to calculate the self-propulsion velocity (|2.9[) 
one needs only to know the main approximation Xi for mutual oscillations of spheres, 
while this approximation is completely described by simple equations (|2.10j) , (|2.11l) . 



2.4. Self-propulsion velocity 

One can see from (|2.12j) that Xq s = I -x s /N = 0(5). Hence the order of magnitude of 
dimensionless physical velocity is: 

Vo = I ■ x t /N = e 2 X Qs = 0{e 2 5) 
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The substitution of A' ()2.2[) and C 1 (|2.13j) into (|2.12p and subsequent algebraic trans- 
formations lead to 

V o = E (2.14) 

i<k<l 

G lkl = 2R i R k R l [ -L + TT - 7T ) felr) (2-15) 

where the sum (|2.14p is taken over all possible triplets (i,k,l) : l<i<k<l<N.ln 
(|2.15|) one can also take 2{l ik l k i T ) = {hk^klr ~ hkrlki), which can be proved by integration 
by parts. The formulae (|2.14[) , (|2.15|) have been obtained for N = 3,4, 5 by explicit 
analytical calculations and for any N by the method of mathematical induction. These 
calculations are straightforward but rather cumbersome to be presented here. However, 
as soon as (|2.14[) and (|2.15j) arc known, they can be verified by separate calculations of 
all terms proportional to XjL 2 ^. for each particular pair i,k. For example, for i = 1, k = 2 
the relevant part of matrix Aq = A' ik (|2.2|) contains only A' 12 = A' 21 = Ri2h/L 2 2 , while 
all other components are zero. The related part of (|2.12j) can be easily calculated; it leads 
to the same expression as the corresponding extraction from (|2.14l) . (|2.15[) . which in this 
case contains only (N — 2) terms with I = 3, 4, . . . , N. 

The number of terms/triplets in (|2.14[) rapidly increases with N: for a three-swimmer 
the sum (|2.14[) contains the only triplet, for a four-swimmer - four triplets, for a five- 
swimmer - 10 triplets, while for a ten-swimmer the number of triplets grows up to 120. In 
general (|2.14p contains N\/[(N — 3)!3!] triplets. It is important to emphasise that (|2.14p 
contains all triplets in a micro-robot, not only triplets of the neighbouring spheres. If 
one takes into account only the triplets of the neighbouring spheres (say, (1,2,3) and 
(2,3,4) out of (1, 2, 3), (1, 2, 4), (1, 3, 4), (2, 3, 4) for N = 4), then they givejDnly the rough 
estimation of Vo> which can be misleading since different correlations (lalpr) can have 
different values (see below). 



Function G123 f|2. 15[) has been introduced by Golestanian & Ajdari (2008) and stud- 
ied by |Alexander et.al. (2009)] |Alouges et.al. (2008)||Golestanian & Ajdari (2009)| in the 



context of a three-sphere micro- robot. We call Giki a Golestanian function. Formulae 
(|2.12p . (|2.14[) . (|2.15[) give the main result of this paper: the self-swimming velocity of an 
A^-sphcrc micro-robot represents a linear combination of Golestanian functions for all 
available triplets. 

2.5. Examples of homogeneous micro-robots: power, velocity, and efficiency 

The explicit formulae f|2.14[) . (|2. 15|) can be used to obtain physically interesting results 
(optimal strokes, required power, related forces, efficiency, etc.) for various TV-sphcre 
swimmers. We briefly address some of these questions below. In all these examples, we 
consider only homogeneous micro-robots, consisting of equal spheres Ri = 1 and of equal 
upstretched arms L a = 1. 

The scalar product of the main equation (|1.6p and x t leads to the average power of a 
micro-robot 

V = (f-x t ) = E 2 (x 2 lT ) + 0(e 2 S) 

where we have taken into account that x t = ex\ T + 0(e 2 ), which follows from (|2.5p 

2 

and (|2.3p . Another expression V s = pV Q represents the power, which is required to 
drag a micro-robot with velocity Vo in the absence of its oscillations (when the main 
approximation for the dimensionlcss Stokes's friction force is — pVo). LighthhTs swimming 



On self-propulsion of N -sphere micro-robot 



efficiency (see Becker et.al. (2003) ) is the ratio £ = V s /V, which in our case is 

P 3 (x 2 lT ) 

where x lT is determined by (|2.f f I) . 

For a three- swimmer , cqns. (|2.14|) . (|2.15p yield: 

V = ^-G 123 , G 12 3^&) (2.17) 

Further simplihcation can be achieved if we accept that oscillations of both arms are 
harmonic and have equal amplitudes 

l a = cos(r + if a ), then 2(hhr) = sin(y>i — ipz) (2-18) 
with constant phases < ip a < 2ir. Substitution of (|2 . 18[) into (|2.17[) and (|2 . 16[) gives 

T7 2c 7 . , c f7eS\ 2 sin 2 

v o = £0—sm(p, E= — — -, <p = ipi-fi 

36 \ 12 / 2 + cos<p 

which shows that their maxima take place at different <p: 

m&xVo ~ 0.19e 2 <5 at = tt/2 ~ 1.57 (2.19) 
max£ = 0.182e 2 <5 2 at = 1.80 

Similar consideration for a four-swimmer yields 

_ s 2 S — — — — 

Vq = — —(Gi23 + Gl2i + G134 + G234) 

16 

which (with the use of (|2.18|) ) leads to 

"^0 = e 2 6— S 4 (<f>, tp), S 4 {<f>, tp) = (l + G)(sin0 + sin^) + 2Csm(<f> + V) 
64 

where C = 41/63 ~ 0.65, <p = <pi — if2, and tp = ip 2 — <f3- Then (|2.16|) can be expressed 
as: 



e = ^ 



7e5\ 2 <S|(0» 



16 / 5 + 2 cos + 2 cos-0 + cos(0 + 
The computations show that: 

maxFo ^ O.Ue 2 6 at ~ V - 1.10 (2.20) 
max£ ~ 0.55e 2 <5 2 at (p ~ tp ~ 1.38 

For a five-swimmer, one can write 

V^O = "^-(Gi23 + G124 + G125 + G134 + G135 + G145 + G234 + G235 + G245 + G345) 

which leads to: 
- e 2 6 

V Q = —S 5l S 5 = (4a + & + c)(sin0 + sinx) + (5a + 26) sin V> + 
50 

(a + 2b + c)[sin(</> + tp) + sin(V> + %)] + (a + 2c) sin(0 + ip + x); 

a = 7/8, 6 = 41/18, c = 151/72; (/) = ^ - tp 2 , tp = tp 2 - <P3, X = ¥3 ~ ¥i 



8 



V. A. Vladimirov 



Then (f2~T6|) takes form: 



1 10) TV 



V5 = 10 + 3 cos 4> + 4 cos ip + 3 cos \ + 
2 cos(4> + ip) + 2 cos(-0 + x) + cos(0 + + x) 



The computations give: 



max Vq 
max£ ^ 



- 0.77e 2 (5 at 
1.00e 2 <5 2 at < 



; x ^ 0.86, V ^ 0.83 
X 1.14, V ^ 1-08 



(2.21) 



(2.22) 



The numerical results (|2.19j> - (|2.22|> show that both maxVo and max£ grow when N 
increases. For reasonably small values of parameters (say, e ~ 0.2 and 8 ~ 0.2) we have 
max£ ~ 0.1%, hence the efficiency of considered micro-robots is low. 

In order to compare the velocities of micro-robots and micro-organisms we use the 
dimensional variables, in which max7 ~ uj*L*e 2 S; it shows that (for typical stroke fre- 
quency of self-swimming microorganisms, which is about several Hz, see |Pedley fc Kessler (1987)| 
Pedley (2009)] |Polin et.al. (2009)[ |Vladimirov et.al. (2004) ) a micro-robot can move it- 
self with the rate about 10% percents of its own size per second. This estimation is 20-^40 



times lower than a similar value for natural micro-swimmers, see Vladimirov et.al. (2004) 



it shows again low efficiency of considered micro-robots. If we suggest that the function 
maxVo(-W) grows with a similar rate (as it has been calculated for N = 3,4,5), then 
micro-robots with TV = 8-^10 could swim with a speed, similar to that of micro-organisms. 



3. Discussion 

1. Our approach (based on the two-timing method and a distinguished limit) is tech- 
nically different from all previous methods employed in the studies of micro-robots. The 
possibility to derive explicit formulae for an iV-sphcre micro-robot shows its strength and 
analytical simplicity. The used version of the two-timing method has been developed in 
Vladimirov (2005 )|| Vladimirov (2008) ~ 



2. One can see that V = Q(e 2 5) (|2T 



Vladimirov (2012b) 



]), which is the same as the result by Golestanian & Ajdari (2008) 



Golestanian & Ajdari (2009) for a three-sphere swimmer. At the same time our choice of 
37 



slow time s = e t 



agrees with classical studies of self-propulsion for low Reynolds 



numbers, see Taylor (1951)| Blake (1971) Childress (1981) as well as geometric studies 



of Shapere fc Wilczek (1989) 



3. The 'triplet' structure of a formula for self-propulsion velocity (|2.14[) can be ex- 
pected without any calculations, on the base of the result for N = 3. Indeed, if we are 
interested in the main term of the order e 2 S, then only triple interactions can be taken 
into account, as they have been described by Golestanian & Ajdari (2008) For example, 
the interactions between four spheres produce a term of the next order 0(e 3 S) in Vq, 
where the multiplier e 3 appears from the motion of three arms involved. 

4. In our examples, all arms move harmonically (|2 . 18|) : it does not provide the maxi- 
mum of Vq. For example, for a three-sphere robot Vo ~ (hhr) (|2.17|) . Since l\ and 1% T 
represent mutually independent 27r-periodic in t functions, it is clear that the maximum 
of this correlation appears when these functions coincide or proportional to each other. 
If max(Zi^2r) is calculated under the constraint of fixed amplitudes (which is natural 
for realistic experimental devices of variable arm lengths), then one can find that the 
theoretical maximum of this correlation is 2/n, which is higher than 1/2 for harmonic 
oscillations (|2.18j) . Such an improvement will increase with the growth of N . In particu- 
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lar, non- harmonic periodic Z a (r), providing the optimal strokes, have been discovered in 



computational studies of four-sphere micro-robots by Alexander et.al. (2009) 

5. In our study we build an asymptotic procedure with two small parameters: e — > 
and S — > 0. Such a setting usually requires the consideration of different asymptotic paths 
on the plane (e, i5) when, say S = 5(e). In our case we can avoid such consideration, since 
small parameters appear (in the main order) as a product e 2 S. 

6. The mathematical justification of the presented results by the estimation of an error 



in the original equation can be performed similar to Vladimirov (2010) Vladimirov (2011) 
One can also derive the higher approximations of Vq, as it has been done by Vladimirov (2010)] 



Vladimirov (2011) for different cases. The higher approximations can be useful for the 
studies of motion with Vq = (say, if all correlations involved in (|2.14[) are zero) . 

7. In the literature quoted in Introduction one can find interesting discussions about 
the physical mechanism of self-propulsion of micro-robots. A clear illustration of this 



mechanism is given by Avron et.al. (2005) At the same time one can notice that the 



self-propulsion of deformable bodies in inviscid fluid represents a classical topic, see e.g. 



Saffman (1967) Kanso & Newton (2009) It is interesting to note that the qualitative 



explanations of self-propulsion in an inviscid fluid and in self-propulsion in creeping flow 
can be seen as the same if one replaces the term virtual mass (for an inviscid fluid) to 
viscous drag (for creeping flows). 



The author is grateful to Profs. A. D. Gilbert, R.Golcstanian, K.I. Ilin, H.K.Moffatt, 
and J. Pitchford for useful discussions. 



REFERENCES 

Alexander, G. P., Pooley, C. M., and Yeomans, J.M. 2009 Hydrodynamics of linked 
sphere model swimmers. J. Phys.: Condens. Matter, 21, 204108. 

Alouges, F., DeSimone, A., and Lefebvre, A. 2008 Optimal strokes for low Reynolds num- 
ber swimmers: an example. J. Nonlinear Sci., 18, 277-302. 

Avron, J.E., Kenneth, O., and Oaknin, D.H. 2005 Pushmepullyou: an efficient micro- 
swimmer. New J. of Physics, 7, 234. 

Becker, D. J., Koelher, C. M., Ryder, and Stone, J.M. 2003 On self-propulsion of micro- 
machimes al low Reynolds number: PurcelPs three-link swimmer. J. Fluid Mech., 490, 
15-35. 

Blake, J. R. 1971 Infinite models for ciliary propulsion. J. Fluid Mech., 49, 209-227. 

Chang, S.T., Paunov, V.N., Petsev, D.N., and Orlin, D.V. 2007 Remotely powered self- 
propelling particles and micropumps based on miniature diodes. Nature Materials, 6, 235- 
240. 

Childress, S. 1981 Mechanics of swimming and flying. Cambridge, CUP. 

Dreyfus , R., Baudry, J., Roper, M.L., Fermigier, M., Stone, H.A. and Bibette, J. 
2005 Microscopic artificial swimmers. Nature, 437, 6, 862-865. 

Earl, D. J., Pooley, C. M., Ryder, J.F., Bredberg, I. and Yeomans, J.M. 2007 Modelling 
microscopic swimmers at low Reynolds number. J. Chem. Phys., 126, 064703. 

Gilbert, A. D., Ogrin, F. Y., Petrov, P.G., and Wimlove, CP. 2010 Theory of ferro- 
magnetic microswimmers. Q.Jl Mech. Appl. Math., 64, 3, 239-263. 

Golestanian, R. AND Ajdari, A. 2008 Analytic results for the three-sphere swimmer at low 
Reynolds number. Phys. Rev. E, 77, 036308. 

Golestanian, R. and Ajdari, A. 2009 Stochastic low Reynolds number swimmers. J. Phys.: 
Condens. Matter, 21, 204104. 

Kanso, E. and Newton, P.K. 2009 Passive locomotion via normal-mode coupling in a sub- 
merged spring-mass system. J. Fluid Mech., 641, 205-215. 

Landau, L.D. and Lifshitz, E.M. 1959 Fluid Mechanics. Oxford, Butterworth-Heinemann. 

Lamb, H. 1932 Hydrodynamics. Sixth edition, Cambridge, CUP. 



10 V.A. Vladimirov 

LAUGA, E. 2011 Life around the scallop theorem. Soft Matter, 7, 3060-3065. 

Leoni, M., Kotar, J., Bassetti, B., Cicuta, P. and Lagomarsino, M.C. 2009 A basic 

swimmer at low Reynolds number. Soft Matter, 5, 472-476. 
Moffatt, H. K. 1996 Dynamique des Fluides, Tome 1, Microhydrodynamics. Ecole Polytech- 

nique, Palaiseau. 

Najafi, A. and Golestanian, R. 2004 Simple swimmer at low Reynolds number: three linked 

spheres. Phys.Rev.E, 69, 062901. 
Pedley, T.J. AND Kessler, J.O. 1987 The orientation of spheroidal microorganisms swimming 

in a flow field. Proc. R. Soc. Lond., B231, 47-70. 
Pedley, T.J. 2009 Biomechanics of aquatic micro-organisms. New trends in fluid mechanics 

research, 1, 1-6. 

Polin, M., Tuval, I., Dresher, K., Gollub, J. P. AND Goldstein, R.E. 2009 Chlamy- 
domonas swims with two gears in a eukaryotic version of run-and-tumble locomotion. Sci- 
ence, 325, 487-490. 

Purcell, E.M. 1977 Life at low Reynolds number. Amer. J. of Phys., 45, 1, 3-11. 
Saffman, P.G. 1967 The self-propulsion of a deformable body in a perfect fluid. J. Fluid Mech., 
28, 2, 385-389. 

Shapere, A. and Wilczek, F. 1989 Efficiencies of self-propulsion at low Reynolds number. J. 

Fluid Mech., 198, 587-599. 
Taylor, G. I. 1951, Analysis of the swimming of microscopic organisms. Proc. R. Soc. Lond., 

A209, 447-461. 

Vladimirov, V.A., Wu, M.S., Pedley, T.J., Denissenko, P.V. and Zakhidova S.G. 2004 
Measurements of cell velocity distributions in populations of motile algae. J. Exp. Biol., 
207, 1203-1216. 

Vladimirov, V.A. 2005 Vibrodynamics of pendulum and submerged solid. J. of Math. Fluid 
Mech. 7, S397-412. 

Vladimirov, V.A. 2008 Viscous flows in a half-space caused by tangential vibrations on its 

boundary. Studies in Appl. Math., 121, 4, 337-367. 
Vladimirov, V.A. 2010 Admixture and drift in oscillating fluid flows. E-print: ArXiv: 

1009.4085vl, (physics,flu-dyn). 
Vladimirov, V.A. 2011 Theory of Non-Degenerate Oscillatory Flows. E-print: ArXiv: 

1110.3633v2, (physics,flu-dyn). 
Vladimirov, V.A., 2012a Self-propulsion velocity of A-sphere micro-robot. ArXiv: 

1206.0890vl, (physics,flu-dyn). 
Vladimirov, V.A. 2012b Magnetohydrodynamic drift equations: from Langmuir circulations 

to magnetohydrodynamic dynamo? J. Fluid Mech. 698, 51-61. 



